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ABSTRACT 


The Risg@ National Laboratory, Roskilde, Denmark, atmos- 
pheric puff dispersion model has been tested for an atmos- 
pheric-marine environment. This three-dimensional model 
Simulates the release of GuaSsian pollutant puffs and pre- 
dicts their concentration as they are diffused and advected 
downwind by a horizontally homogeneous, time-dependent wind. 
Atmospheric characteristics such as turbulence intensity, 
potential temperature gradient, buoyant heat flux and maxi- 
Mum mixing depth have been considered. Model predicted pol- 
lutant concentrations have been compared to airborne sampled 
observations. The effect of coastal turbulence not observed 
by the single point meteorological measur2ments made onboard 
ship greatly affects the advection and diffusion of a plume 
as it moves onshore. Additional measurements/predictions 
particular to the coastal area will have to be incorporated 
into the model for it to accurately predict the onshore 


movement of pollutants. 
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I. INTRODUCTION 


The downwind transport and distribution of atmospheric 
pollutants from an isolated source over land or water has 
become an important environmental factor in today's society. 
The need to understand the distribution of smoke, unpleasant 
or potentially harmful foreign gases and perhaps radioactive 
debris from a nuclear powerplant accident are becoming more 
and more essential for industral operations and construction 
planning. The dispersion of such atmospheric pollutants is 
commonly modeled by a standard Guassian plume model which 
computes one-hour average plume characteristics. 

The Meteorology Section of the Ris¢g National Laboratory, 
Roskilde, Denmark, has recently developed a puff model for 
prediction and simulation of atmospheric pollutant diffu- 
S20n. The model considers individual puffs of pollutants 
with specific release rates that are advected by a horizon- 
tally homogeneous wind over a grid. The wind input may be 
either the measured wind from a singl2 point, a spatial 
average ora wind simulation. The model simulates the 
instantaneous plume characteristics by adding ae group of 
puffs, growing in size, as they advect with the wind. A 


Guassian plume nodel, on the other hand, provides a time 


10 





averaged concentration pattern based on a single time aver- 
age wind vector. In the puff model, the plume advects with 
a time series of actual wind data. Thus, the puff model is 
able to predict time varying concentration distributions in 
actual changing wind conditions, making it an appropriate 
tool for dynamical computations of downwind dispersions of 
pollutants. 

A basic comparison of a puff model simulation and a typ- 
ical plume is illustrated in Pig. 1. Looking from above, 
the instantaneous behavior of a plume being advected from a 
source by the wind is shown. The outer cone-shaped contours 
represent the outer limit of the plume boundary and are 
identical in both Figs. 1 (A) and (B). 

Ed 1(A) shows an instantaneous depiction of an actual 
plume. The long-term average plume concentration is shown 
on the extreme right as a smooth curve with a maximum on the 
central axis. Also shown is the instantaneous volume concen- 
tration considered realistic but is of such a_ short tine 
scale that it cannot be predicted or easily measured. 

The puff model prediction is depicted in Pig. 1(B). The 
Circles show the boundaries of individual puffs of pollu- 


tants released from “the source. These puffs are advected 
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Figure 1. Instantaneous Behavior of a ipa Plume and a 
Series of Puffs from a Puff Model 


and diffused downwind by a frequently updated wind. The 
long term average concentration prediction of the puff model 
is expected to be identical to the long term concentration 
of Fig. 1 (A). The short term average pollutant prediction, 
@ Gaussian curve shown on the extreme right, 1s not com- 
pletely realistic but is a reasonable approximation to the 


instantaneous plume concentration profile. 
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The purpose of this thesis is to evaluate and test the 
general characteristics ‘and capabilities of the Ris¢ puff 
model. The adjustable input data will be varied and pre- 
dicted results for various input combinations compared. A 
preliminary comparison will also be made between predicted 
results and data collected from a coastal region using 


observed meteorological forcing data as model input. 
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II. RIS¢ PUFF MODEL 


A. GENERAL CHARACTERISTICS 

The Ris¢g Puff Model is a three-dimensional computer 
model used for the prediction and/or simulation of the dif- 
fusion and advection of atmospheric pollutants. The puff 
model technique is to simulate a plume with Guassian shaped 
puffs with specified release rates within a specified grid. 
The initial size of the puffs is normally one meter in dian- 
eter although this can be easily adjusted. The amount of 
material in a puff is the release rate times the elapsed 
time between puffs. Therefore, a long elapsed time between 
puff releases results in a higher initial puff pollutant 
concentration than a short time interval. This should not 
normally be of concern if an adequate balance is maintained 
between grid size, advection speed and puff release rate. 

The location of the puffs on the grid is determined by 
computing their movement for a finite time step using a 
measured wind field. The growth and buoyancy of the puffs 
are computed from simultaneous specifications of atmospheric 


turbulence intensity and stability and from buoyant heat 


14 





flux at the source. An inversion cap through which pollu- 
tants cannot pass and the source height where pollutants are 
released are variable and can easily be adjusted. Grid dis- 
tances within the model may vary from meters +o kilometers 
and time durations from seconds to hours are possible. 

This puff model has the capability of monitoring a maxi- 
mum cf twenty-five sources of puffs and its grid may contain 
up to 100 puffs. A puff source can be Located anywhere on 
the grid and have a unique relsase rate, start and stop of 
release time, and heat production. When the center of a 
puff moves outside the boundaries of the grid (either hori- 
zontally or vertically), that particular puff is dropped 
from memory. In this way the model does not store irrele- 
vant puff information, thus ke2ping computer memory require- 
ments to a minimun. 

A variable to control the amount of reflection/absorp- 
tion of the pollutant by the surface is easily adjusted in 
the puff model. Such a parameter is of great value both in 
actual dispersion problems and also for gaining understand- 
ings of the plume/surface relationship. 

The model calculates the concentration at each grid 


point by summing the contributions from surrounding puffs 
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for each advection step. The grid concentrations can 
be allowed to accumulate or simply be updated with the lat- 
est instantaneous value. A minimum grii concentration of 
interest can be set to reduce computer run time by dropping 
concentrations too small to be of interest. 

The output of the model contains psriodic results of 
puff locations and concentrations as well as initial input 
verification. The time interval for the periodic results is 
adjusted by the input data. This recurrent lineprinter out- 


put contains: 


e X-Y plane plots showing the position 
of the sources and of puffs inside the 
re f.an! [alr 

e X-Z plane plots of puff positions 
for evaluating plume rise for 2ach 
vertical level of interest, and 

e a table listing of the grid point 


concentrations for each level. 


A computer drawn contour chart of the magnitudes of the pol- 


lutant concentrations is also available. 
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When considering distance between gridpoints (delta 
Moy, 2), only spatial resolution and computer resources need 
be considered. Calculated concentration accuracy is not 
related to the grid-point separation. fo insure that no 
essential information on individual puffs is "hidden" 
between grid points, the grid separation should be adjusted 
dependent upon the size of the puffs at the downwind dis- 
tance of interest. Other specific model configuration con- 
siderations are described in the following sections. They 
are also discussed in more detail in the model behavior 


chapter. 


B. WIND FIELD 

Once a puff is released, it is advected based upon wind 
data measurements at a single point only, normally the 
release point. This limits the validity of the model to 
Situations where the wind field and turbulence can be 
assumed to be horizonally homogeneous throughout the grid. 
It is therefore important to insure that the data obtained 
from such a single point measurement is representative of 
the wind structure for the whole area of interest. 

The wind data are normally obtained in the form of a 


horizonal velocity time series. A vector sequence is formed 
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by averaging over a convenient interval. These data are 
read into the model after being segregatei into turbulence 


classes as discussed in the next section. 


Ce TURBULENCE INTENSITY AND DIFFOSION 

The growth/diffusion of a puff depends upon the turbu- 
lence intensity. To account for this growth, the puff model 
applies the theory for relative diffusion suggested by Smith 
and Hay (1961). 

The turbulence intensity is defined to be the standard 
deviation (Sigma) of the wind direction (in radians) 
squared. The sigma values are collected for the same short 
time periods as the Se ucda measurements used to advect 
the puffs. Therefore, the intensity of the turbulence which 
governs the relative diffusion of the puffs, can be adjusted 
along with the the advecting wind speed after each time 
step, if conditions warrent. 

A very low value of turbulence intensity (as 0.0002) 
represents a small standard deviation (0.9), normaily a sta- 
ble atmosphere and a weak puff dispersion/diffusion. As the 
atmosphere becomes more unstable, the turbulence intensity 
increases along with an increase of sigma values and plume 


dispersion sdiffusion. While these characteristics are 
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resentative of turbulence over land, they can be applied to 


over water cases in a broad sense. 


De PLUME RISE 

In the vertical direction, puff-ris? can be accounted 
for by Briggs (1970) plume rise theory. In this case buoy- 
ancy is assumed to be conserved (adiabatic motion), and 
pressure forces, molecular viscosity and local density 
changes are considered small and are neglected. The rate at 
which a puff rises as it is advected downwind is a function 
of the buoyancy flux, wind speed, puff distance traveled and 
stability of the atmosphere. Plume rise is considered sep- 


erately for each individual puff. 


Ee REFLECTION 

The interaction of the pollutant with the surface is 
adjustable and can be easily changed in the input data. 
Total reflection or absorption or a fraction between the two 


can be used. 


re LIMIT OF MIXING DEPTH 
The effect of an atmospheric lid (inversion) can be 
applied in the model to limit the vertical movement of the 


pollutant. The model does not permit the plume ‘+o rise 


19 





above this cap. When a maximum mixing level is in effect, 
it acts to totally reflect the pollutants in the same manner 
as total reflection at the surface. This mixing cap would 
act aS an inversion when the puff would be expected to grow 
much more readily in the horizonal than in the vertical 


direction. 
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III. DATA COLLECTION 


An intensive field tracer study was p2rformed during the 
fall of 1980 and winter of 1981 in the Santa Barbara Channel 
area of the California coast. The work was supported by the 
Bureau of Land Management (BLM) and performed by the Envior- 
mental Physics Group of the Naval Postgraduate School (NPS) 
and AeroVironment Inc (AV), Los Angeles, CA. This study was 
designed to help validate and/or modify Saussian dispersion 
models for coastal use and to build a data base for future 
model development. Air pollution models in current use have 
not been adequately validated for the over-water regime. 

In the experiment, a tracer gas (SF,) was released sev- 
eral miles offshore from the NPS Research Vessel (R/V) Aca- 
nia. Ambient gas concentrations as low as 10 parts per 
Stallion (PPT) were determined by an array of land based 
sensors, from a small boat and at various levels by an air- 
craft equipped with a continuous SP, analyser. A chart of 
the experimental area and locations of th2 various platforms 


momshown in Fig. 2. 
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Figure 2. 





The aircraft flew through the plume at various eleva- 
tions offshore and overland. The plume transect tracks per- 
tinent to this study were made parallel ts the coast approx- 
imately one-half mile offshore. The airborne sampling 
consisted of instaneous concentrations (PPT) at selected 
points at different levels over a period »2f six hours. The 
observations, recorded at locations 10-70 at altitudes of 61 
and 91m above MSL on January 29, 1980 ar2= shown in Fig. 3. 
Average concentrations over the noted tine period are shown 
at the bottom of each altitude block. 

The following marine meteorological parameters were 
measured onboard the R/V Acania while anchored approximately 


7.4 km offshore: 


erelative wind speed eair temperature 
ewind speed fluctuation edéew point 

esea surface temperature eShip roll 

esky cloud cover eship location 


erelative wind direction 
einversion height (acoustic sounder) 
evertical temperature and humidity profilies 


(Shipboard radiosonde launch every 12 hours). 
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DATE: 9/29/80 TIME: 1700-1800 POT 
TRANSECT ALTITUDE: 300 FEET MSL 1/2 MILE OFFSHORE 


BEG END 
17:20:50— 17:21:14 


-10 288 355 488 56 
-142 521 288 421 156 


305 471 1634 372 


17:33:35- 17:34:04 
471 903 1169 39 
185 628 219 82 20 
374 506 174 -26 
30 | 40 
DISTANCE (GRID COORDINATES) 


17:20:50- 17:59:42 





DATE: 9/29/80 TIME: 1700-1800 PDT 
TRANSECT ALTITUDE: 200 FEET MSL 1/2 MILE OFFSHORE 


BEG END 
17:25:43-17:26:18 





1342 4967 4967 1753 #707 595 184 
7O7) 3024 5304 3921 520 558 109 


2129 7249 6053 3848 1419 372 166 
=1-5 335 5081 7922 2465 634 634 335 


17:37:00— 17:37:41 


671 3548 6108 3903 2277 1007 276 93 


17:25:43- 17:37:41 
316 1679 5192 5821 260 1511 #1371 317 167 


30 40 
DISTANCE (GRID COORDINATES) 


Figure 3. 
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The tracer gas was released at a fixed rate through the 
exhaust of one of the ships motor generator sets. The gen- 
erator was run at a constant speed resulting in a constant 
stack temperature and flow rate. 

The above collected meteorological information enabled 
the data grid to be established and th atmospheric wind 
field, turbulence intensity, stability and buoyancy flux to 
be derived for inputs into the puff model. The puff model 
prediction based on this actual data formed the basis for 
the model performance evaluations carried out here. From 
this basis, different input variables were adjusted to note 
the effect on the advected concentrations--both in relation 


to each other and to the airborne measurements. 


25 





- IV. MODEL BEHAVIOR 


The Puff Model was run on the NPS IBM S/370 Model 3033 
AP computer with two goals in mind; (1) familiarization with 
model performance wnder actual conditions and (2) a compari- 
son cf model predictions with observed data. The two goals 
are interrelated in the sense that atmospheric data col- 
lected in the aforementioned tracer study were used to form 
an initial prediction of the plume dispersion, and varia- 
tions of that data were used to evaluate the limits of the 
model. The data used as. input to the model represented the 
Marine atmospheric conditions as determined from R/V Acania 
meteorological data at the time of the exoeriment. 

Proper grid spacing was arrived at by considering puff 
spread, mean wind direction and the geographical area of 
interest. With the initial prediction in hand, data input 
variables of the model were adjusted and their effects 
(changes in prediction) noted. All model predictions were 
compared with the aircraft observed data. 

A724 X 4.3 km downwind area of intarest was initially 
gridded into an 17 X 10 array. Distances between horizonal 


and vertical grid points were approximately 435 no (Pig. ZN 
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Since many of the aircraft observation times centered 
around 1730 hours (all times are Pacific Daylight Time), 
model puff releases were initiated at 1630. The 30-minute 
wind speed averages obtained from data taken onboard the 
ship between 1630 and 1730 were 4.7 and 4.8 ms~4t. The first 
advected puffs would be expected to arrive at the back edge 
of the grid slightly before 1700, and by 1730 a steady con- 
Sistent plume would be passing through the area of the air- 
craft track. (The model showed, in fact, puffs leaving the 
back edge of the grid slightly before 30 minutes after puff 
release). 


The average wind direction was recorded on the ship 


every 15 seconds. The standard deviations of the wind 
direction (sigma) were computed as approximate one minute 
averages. These, in turn, were averaged over 30 minutes to 
correspond with «he 30-minute wind speed averages. The 


Sigma values during the time of interest were 1.0 and 0.9 
resulting in the very small turbulence intensity values of 
-0003 and .0002. 

A delta Z value of 33 m was used to observe plume con- 
centrations at the altitudes of 0, 33, 66 and 99 m above the 
surface. These levels were chosed for tomparison with the 


aircraft transect altitudes of 61 and 91 n. 
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Fine scale vertical temperature and humidity plots were 
drawn based upon radiosonde soundings taken onboard the 
ship. The sounding taken at 1735 PDT (Figure 4) shows a 
shallow unstable layer near the surface topped by an inver- 
sion extending to near 400 mn. A 80 m depth of the mixing 
layer was subjectively established. The potential tempera- 


ture gradient computed by the formula 


a aor 
s- = <> + .0098 z (1) 


Qo 


was found to be 1.0 deg K/100 nm. 

Basic data to determine source strength and heat emis- 
Sion from the ships stack were taken from Schacher, et al 
(1981). As previously mentioned, SP, gas was released 
through the ships motor generator exhaust ata constant 
rate. The stack temperature was 250 deg F, the flow rate 
was 7.13 *® 103 m3s-t. The SF, release rate was 47.91 1b 
Pe !. The top of the ships stack--considered +o be the 
source 2levation-- was 4 mn. The source strength was con- 
verted to 6.04 gm s~! for input into the nodel. 


Heat emission (H) in KW was determined by the formula 
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et a) «(SSCP OW (2) 


where 


delta T = temperature (stack - air) d3g K 


R = dry air gas constant 
= 2.87106 erg/g deg K 
P = 103 * P(mb) = dyne/cmn? 
Cp = specific heat of dry air 


-24 cal/g deg K 


Flow = 16.39 *® 7.13 * 103 cn3/s 


The heat emission was thus computed as 15.07 KW. 

Initially the model grid was established after noting 
the area of maximum airborne sampled concentrations (between 
points 24 and 43 of Pig. 3) and the location of the ship. 
It became obvious during early model runs that, with the 
actual wind direction input, the model predicted plume was 
being advected south of the grid towards point 60 on the 
aircraft track. Obviously, the steering wind, as measured 
onboard ship, was not constant 211 the way to the shore. A 
northward turning of the plume was indeed detected several 


times during the experiments by the aircraft. To compensate 
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for this effect, the source of the plume release was moved 
in the model three grid spaces (1305 m) to the north so that 
the maximum predicted plume concentrations would pass 
through the areas of the maximum airborne measured concen- 
trations. No corrections wer2 made to the model predicted 
plume concentrations because of this adjustment. However, 
one could reason that the predicted concentration values 
would be higher in comparison with measur2d values since the 
coastal turbulence and wind shift--which would tend to dif- 
fuse the plume--were not considered. 

Initial model runs with the small turbulence intensity 
classifications of .0003 and 0002 failed to show plume con- 
centrations greater than 1 ® 10-!2 iain th2 grid at any level 
Other than at the source. Apparently, th2 grid spacing was 
too large andthe narrow plume was advecting between the 
grid points. In an effort to locate the plume, a combina- 
tion of model runs were performed varying the turbulence 
intensity and grid spacing as shown in Table I. 

In this table, a mixing level cap of 80 m was in effect 
for the model predictions. No concentrations above that 
level were allowed in the computations. As previously men- 
tioned, the model mixing level cap totally reflects all 


pollutants back downward. 
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With a grid spacing of 435 m anda turbulence intensity 
of .05, a plume concentration covering four grid spaces at 
the east end of the grid was produced. Predicted plume con- 
centrations slowly decreased as the turbulence intensity was 
reduced to .02 (atmospheric stability increased). At an 
intensity of .01, the concentration dropped by about five 
orders of magnitude. Normally, one would expect increased 
concentrations with increased atmospheric stability. Per- 
haps the result noted here is due to the plume shrinking 
away from a grid point (and becoming more concentrated 
between the recorded grids) with the increase in stability. 

To increase the grid resolution, the grid spacing was 
reduced by half to 217.5 m and again to 108.75 m. With each 
reduction the grid was reduced by half in the "Y" direction 
and doubled in the "X" direction thus keeping distances 
between grid spaces equal in all directions. This of course 
greatly increases the computational requirements. Tiron Ly 
plume predictions along the back edge are needed (as in 
Table I), the downwind grid distance may be held constant at 
435 m while the horizonal crosswind resolution is increased. 
In this way many unnecessary computations are not made. 


However, the increased horizonal crosswind resolution is 
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computed over the entire downwind grid, which in this case, 
is not necessary. A more satisfying solution to this prob- 
lem is to install the capability of using a variable resolu- 
tion grid with the model so that downwind areas of 
particular interest can be covered with 2 dense grid while 
other areas of not so much interest can be sparsely grided. 
In order for the advection of the plume to remain on the 
array when increasing the horizonal crosswind resolution and 
decreasing the area exposed on the grid, the plume source 
was adjusted along the western boundary of the grid. The 
relationship of the vertical grid points to changes of the 
source location is shown in Table II. The plume source for 
each grid resolution is noted with an arrowhead. Grid 
points that are aligned vertically in the table have identi- 
cal locations and should have the same predicted plume con- 
centrations. As mentioned earlier, an increase of grid 
resolution does not affect the predicted concentration. 
Notice that for the same grid points in Table I, the pre- 
dicted concentrations with a turbulence intensity of .01 
remain constant with changes in grid spacing--only the grid 


resolution was changed. 
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TABLE IT 


The Relationship of the Y Axis along the Western Grid Edge 
. to changes of Grid Distance between 435 and 108.75 nm. 


DELTA 


Y (m) Grid Points Along the Y Axis 

435 4 5 6 7 8 9 
217.5 0 1 Z 3 4 5 6 7 8 9 
108.75 O 123456 7 8 9 


A 


From Table I it is obvious that the 435 m grid spacing is 
too large and that the higher resolution does indeed "seet 
concentrations that would otherwise be missed. 

The problem of a increasingly narrow distance covered on 
the grid as resolution 1S increased can sometimes be at 
least partly corrected by reversing the K and Y coordinates 
and adjusting or rotating th2 advecting wind direction. 
This can easily be done with the use of the "TURN" model 
input parameter. This procedure sometimes becomes necessary 
Since one of the grid directions is limited by the width of 
the output printer paper to less than or equal 10 grid 
oats. 

The model input variables, meteorological and source 
values were adjusted to note their effect oon plume concen- 
trations. A deeper understanding of how the model works and 


how the atmosphere affects dispersion can also be gained by 
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such adjustments. A turbulence intensity class of .05 was 
used, except when studying intensity itself, because it had 
previously demonstrated a good iownwind grid coverage of the 
plume. 

In order to note the effect of the maximum mixing level 
on the plume concentrations, saveral model predictions were 
run, varying only the height to which the plume was allowed 
to rise. Exact grid point reproductions were not possible 
Since the model only allows the height of the mixing level 
to be an integer multiple of delta Z. The vertical grid 
Spacing is therefore not equal. However, the antipicated 
trend of increased concentrations as th2 mixing level is 
lowered is evident from Table III. 

The reflection/absorption of the plum2 at the surface is 
controlled by the model variable "REFLEC". Tests of the 
extremes oof total absorption (0.0) and total reflection 
(1.0) were performed. The results showed a 50 percent 
reduction in plume concentrations at the west end of the 
grid with total absorption compared to total reflection in 
otherwise identical model runs. 

The model has a self-imposed limitation cf 100 puffs 


from all sources on the grid. The model will terminate if 
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TABLE I1fifI 


Plume Cone-Prpatlors between Surface and 99 m under 
Different Maximum Mixing Levels. 


Max Mixing 
Level (5) (6) (7) (8) 
None -44E-8 -55E-5 ~20E-5 ~ 225-9 Stake. 
-43E~8 -54E-5 ~LIE=-5 ~21E-9 33 m 
-40E-8 -50E-5 -L18E=-5 -20E-9 66 m 
- 36E-8 ~45E-5 - 16E-5 -18E-9 99 m 
80 m -57E-8 ~72E-5 ~20E-5 -28E-9 Sfc 
-59E-8 ~ 74E-5 ~27E-5 -29E-9 40 m 
-58E-8 ~72E-5 -26E-5 -29E-9 80 m 
0 0 0 0 > 80 m 
30m -64E-8 ~8lE-5 ~29E-5 -32E-9 Sfc 
-65E-8 -8lE=-5 ~29E-5 ~32E-9 30 m 
0 0 0 0 >30 m 


this number is exceeded. A balance must be made between the 
rate at which puffs are released from th2 source (TAU) and 
the time it takes the puffs to be advected across the grid. 
A release rate of one puff every 40 seconds was predomi- 
nantly used during this study. 

The turbulence intensity variable was varied to include 
conditions that are more unstable. As atnospheric instabil- 
ity increaseés, the plume would be expected to expand 
whereas, with stable conditions, the plume should remain 


narrow and highly concentrated. 
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With the use of NPS contouring routines and the 
subroutine “DRAG, a visual comparison of the plume disposi- 
tion and concentrations is available. Since the computed 
plume concentration varies ovar many orders of magnitude, 
the concentration values were converted to integer numbers 
by multiplying by 1.2 ® 10!3 and then taking the logarithn. 
These logarithms are then smoothed. EAWwS, a contour plot 
representing order of magnitud2 concentrations was be pro- 
duced. As with the model variables "MAPTIM" and "KPLANS" 
which control the frequency and vertical levels of printer 
plots, "DRAW" can be called to contour concentrations at any 
time period and for any level required. 

Plume concentration distributions for turbulence inten- 
Sities of .05, .10 and .25, all other variables cons<ant, 
are shown in Figs. 5, 6 and/7. As expected, the plume 
becomes wider and the concentration decreases as the <urbu- 
lence intensity/diffusion increases and the atmosphere 
becomes less stable. 

In Figs. 5-7, the plume source was located at grid point 
"0 6). The wide plume in the lower part of the plot is not 
real but is a function of th2 smoothing routine spreading 


the early puffs more than would be expacted. Since the 
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Figure 5. Orders of Magnitude of Plume Congentration with 
Turbulence Intensity equal to .05. 
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Figure 6. Same_as Figure 5 except [Turbulence Intensity 
equal to .10. 
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Figure 7. Same_as ice > except Turbulence Intensity 
equal to .25. 
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smoothing routine 
near the source, 
primary interest is 


expect the puffs to 


(1979), showing the 
concentration, the 


advecting wind speed 


PUFF SIZE 


the smoothing 


would tend to smooth strong concentrations 


should be eliminated if the 


near the source. Actually one would 


behave as in Fig. 3, from Mikkelsen 


relationship between the puff size and 


rate of puff release (TAU) and the 


U. 





CONCENTRATION 





Bgure 8. 


Relationship between Puff Size, C 
Puff Release Rate (TAU) and Advecting Wind Speed 
GO (Mikkelsen, 1979). 


oncentration, 
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Puffs would have to travel the distanc? Xmin before they 
expand to a size where they effectively overlap and form a 
solid plume. From Figs. 5-7, one can see that plume concen- 
trations have increased with distance and that a Xmin has 
been reached in the middle to upper part of the grid. It is 
at this point that the puff mod21 would be expected to accu- 
rately predict plume concentrations. If the area of inter- 
est is before the present Xmin, the release rate of puffs 
would need to be increased so that th successive puffs 
would overlap sooner. Also noted that as the turbulence 
intensity increases, the area of maximum concentration of 
the plume expands while the central concentration decreases. 
This agrees with conservation of mass theorey. 

To appreciate the relativ2 importance of the source 
strength and buoyant heat flux, these variables (discussed 
in Chapter IV) were doubled separately and together and the 
concentrations compared to the concentrations from the 
actual conditions. Little or no changes in concentration 
were noted when the buoyant heat flux was doubled and source 
strength remained the same. However, when the source 
strength was doubled and heat flux held sonstant, the grid 


concentrations doubled as expected. Thus, under existing 
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conditions, the source strength was critical to the pre- 
dicted plume concentrations while the buoyant heat flux, 
within the range tested, was not relevant. The vertical 
printer plots did show an initial puff rise soon after 
release due to the initial heat release but as the puff rose 
and expanded, it soon reached the ambient temperature and 
leveled off. The buoyant heat flux would probably be more 
important when dealing with a smaller scale grid or greater 


heat release. 
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No attempt was made to compare actual puff model concen- 
+ration predictions at exact grid points to aircraft obser- 
vations for the following reasons: 

e The aircraft locations were approximations--the exact 

locations were not known. Large differences in predicted 
concentrations are seen with small grid separations as 


evidenced in Table II. 


e AS noted in Fig. 2, the aircraft observations were taken 
over a period of time at different levels--while the puff 


model produced multilevel instantaneous predictions. 


e AS mentioned earlier, the actual wind was not constant 
between the ship observation site and the oppsite side 
of the grid near shore. Since the model advects the 
puffs based upon ship observed wind, the behavior of 


actual plume would be different from predicted. 


e Calabration procedures for the SFg continuous analyzer 
mounted onboard the aircraft were not available for 


instantaneous concentrations greater than 1010 PPT. 
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Therefore, a question of actual levels of SF, 
concentration in the higher ranges exists. 

The puff model predicted concentrations are expressed in 
g/m3 while the aircraft observations are shown as the volume 
of SF, per wnit volume of air in PPT. A conversion between 
the predicted and observed concentrations was obtained by 
computing the partial pressure and molecular weight of SF, 
at standard pressure and temperatur2. A conversion of 

g/n3 = (.63%10-11) & observed concentration (PPT) 
was thus found. 

Generally, the aircraft sampled concentrations (Fig. 2) 
show values between 100 and 8000 PPT. Converting these 
observed concentrations to predicted concentration units 
gives values between .63%10-9 and .50%10-7 g/m3. These 
observed concentrations are much smaller than the values 
shown in Table II. Perhaps this difference could be 
explained by the fact that the puff model advected the plume 
toward the coast in the same direction under the same very 
stable conditions as observed on the ship. Any considera- 
tion of increased turbulence and wind shifts near shore 
would be expected to reduce the actual plume concentrations 


toward the observed concentration levels. 


46 





Increasing the turbulence intensity to 0.25 and keeping 
all cther variables constant, the concentration values would 
decrease to the order of magnitude of 1)-®--closer to the 
observed concentrations. (This would require the wind 
direction standard deviation to increase from 1to 28). 
However the increased instability would cause the plume to 
spread over a much greater area (Fig. 7) than observed by 


the aircraft. 
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VI. CONCLUSION 


AND RECOMMENDATIONS 


The puff model has been demonstrated to be a versatile 
working dispersion model. Different combinations of input 
variables showed the expected reasonable results. The dif- 
ferences between model predicted and aircraft observed plume 
concentrations do not seem to be the fault of the model but 
mainly that of the highly variable meteorological conditions 
found along a coast. 

Probably the most obvious conclusion reached from this 
study is that predicting the behavior of 2 plume moving over 
@ marine environment onto a coastal regidn has significant 
problems. In all probability, atmospheric boundary layer 
conditions 7.4 km offshore can be very different from those 
observed in the more turbulent coastal region. The single 
point meteorological measurement at the sdurce2 should not be 
expected to adequately represent plume characteristics as it 
nears a meteorologically variable coastline. Addie ional 
observations (primarily wind sp2ed and direction), or other 
means of predicting the coastal meteorological conditions, 
would have +to be incorporated into the puff model to ade- 


guately handle this problen. 
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The advantage of incorporating variable grid spacing 
within the puff model and the obvious benifits have already 
been discussed. 

Presently, the mixing cap of the puff model is required 
to be located at an integer multiple of dalta Z. More flex- 
lability in this parameter to include any level, regardless 
of delta Z, would be benifical. 

Along with the puff locations shown on the lineprinter 
output, a maximum concentration level of ¢ach puff would be 
helpful. 

In future experiments, several aircraft tracks should be 
made further out from the coast in an attempt to avoid the 
turbulent coastal region. Obs2rvations thus obtained in a 
noncoastal environment would help to varify the model prs- 


dictions without the coastal influence. 
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APPENDIX A 


MAJOR SECTIONS OF THE PUFF 4ODEL 


The Risg Puff Model has been described by (Mikkelsen, 


1979). The code also is well documented with comment state- 
ments. With that information and the outline to be provided 
in this and the following appendices, the computational and 
input/foutput procedures will be obvious. 

The program and input data are stored on cards for the 
sake of permanency. For efficient operational execution, 
the program and input data cards are read on a disk within 
the computer. The model can then be run at will without 
reference to the original data cards. Minor changes can 
easily be made directly on the disk both to the model and/or 
data before each execution. 

The model can be separated into the following main sec- 
tions: = 

A. Input Data 

B. Initial 

Ge Calculating 

D. Output 

E. Error Diagnostics 


F. Subroutines 
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These will be described separately in the following sec- 


tions. 


Ae INPUT DATA SECTION 
The input data includes the variables shown in Table IV. 


TABLE IV 
Input Data Variables for the Puff Dispersion Model. 


Wind History Potential Temperature Gradient 
Turbulence Intensity Buoyant Heat Flux 

Grid Dimensions Minimum Concentration of Interest 
Mixing Depth Reflection at Ground Level 


Source Locations, Start/Stop Time, Strength, Heat Emission 


Number of Seconds between Advection Steps 
Number of Seconds between Printouts/Plots 
Number of Seconds between Puff Releases 


' The wind field and stability class for th current time step 


are read at the start of the calculation section. 

The variables listed above are printed as a input data 
check and a permanent record to accompany the actual output. 
In most cases the print command can be overridden by YES/NO 


options. 


B. INITIAL SECTION 
Based upon the input data from section (A), the initial 


section specifies and initializes parameters tc be used in 
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the calculating section and is passed only once during =xe- 
cution of the model. The grid and some counters are ini- 
tialized. Constants relating to reflectance, mixing depth 
and stability as well as those controlling the size of some 
of the loops within the model are established. Parameters 
such as number of puff releases per second, number of advec- 
tion steps per second and number of advection steps per puff 


release are determined. 


C. CALCULATION SECTION 

Using current wind and stability class data read at the 
start of the calculation section, the model advects the puff 
centers and calculates ange growth rate and plume rise of the 
puffs. It removes the puffs that have left the grid (hori- 
zonally and/or vertically). The predicted concentration is 
computed at the grid points to include pollutants from all 


rearby puffs. 


D. OUTPUT SECTION 
For time intervals designated by the input data, printer 
plots of the xX-Y¥ and Y-Z grid are producei. A maximum mix- 


ing level is marked on the Y-Z grid if in effect. 
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These plots include the source location and a trace of the 


plume from the release time to the maptime. Also printed at 
this interval is a X-Y table of grid concentrations for each 
vertical level of interest. These concentrations can be 
either accumulated or actual concentrations at the plot 
time. 

Added to the puff model is a versatic plotter routine to 
smooth and contour the grid magnitude concentrations of the 


above tables. 


E. ERROR DIAGNOSTIC SECTION 

If the model is directed by the inpu+ data beyond the 
limits of the design of the program, the program is termi- 
nated by way of the error diagnostic section. Ee. PpELots 
comments relating to commonly made input errors enabling the 


user to isolate problems. 


F. SUBROUTINES 

The subroutine "Sigris" calculates th2 puff size in the 
horizonal and vertical directions. It also estimates plume 
rise associated with pollutent buoyancy. 

The subroutines "Ispace" and "Rspace" are used in the 


framework of the printer plots. 
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The subroutine "Draw" 


to a logarithmic values, 


converts the plume concentrations 


smoothes and 


using NPS inhouse contour subroutines. 


verted to their logaritha 


Values so that 


touring over many orders of magnitude is 
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then contours then 
The values are con- 
the problem of con- 


simplified. 





APPENDIX B 


PUFF MODEL FLOW CHART 


INPUT 
SECTION 





READ - data input variables (primary, source, 
intensity, current stability); averaging 
' timestep (tav) for wind data 


INITIAL SECTION [initialize constants | 


READ & STORE all wind speed and direction 
for current intensity class 


CALCULATION | next wind vector 
RECALL from storage one wind vector for next tav 


@ ae C4 e = OT sten 
Release puff, execute advection step, 
calculate concentration 


> advection 


SECTION See. 
eG. Vv 












. OUTPUT 


SECTION 
PUL pest eters, 


+ yes grid concentra- 


| tions, contour 
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APPENDIX C 


PUFF MODEL CONTRACTIONS 


CHEMIN--Minimum grid concentration of interest 

DELX ,DELY, DELZ-~Distance in meters between grid points 
DOSE--Allows the concentration natrix to accumulate 
DTDZ--Potential temperature gradient (K/M) (.GE. 0) 
HEAT--Individual source heat emission (KW) 
ICOLS--Number of columns in grid (.<LE. 10) 
INST--Instantaneous concentration matrix 

ITIME--Start time 

JROWS--Number of rows in grid 

KPLANS-~Number of vertical levels in gri includes surface) 
MAPTIM--Number of seconds between printer plots 
NRELSE--Number of seconds to stop of release 
NRMULT--Number of sources (.LE. 25) 


NTADV--Integer number of seconds between advection steps 
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REFLEC--Reflection at ground level (0. none;1.0 total) 
SOURNR-- Number to identify source 

SOURST--Strength for individual source (gn/s) 
STOPRL--Individual source stop time (s) 
STRTRL--Individual source start time (s) 

TAU--Integer number of seconds between puff releases 
TURN--Angle of rotation of wind direction 

XSOURCE--X coordinate of source in grid units 
YSOURCE--Y coordinate of Bees in @gead- units 


ZM--Limited mixing depth (n) 
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